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I.  Introduction 

A  plasma  with  hypothetically  infinite  size  is  a  kind  of  media  for 
electromagnetic  waves,  and  their  dispersion  relations  in  various  plasmas  have  been 
extensively  studied.1'3  Several  schemes  of  plasma  production  and  heating  using  the 
waves  have  been  proposed  and  investigated  experimentally  and  numerically;  plasmas 
for  fusion  experiments  as  well  as  processing  plasmas  are  controlled  using  launching 
antenna  and  high-power  wave  sources.  One  of  the  fundamental  facts  we  have  learned  in 
microwave  plasmas  is  that  microwaves  cannot  usually  propagate  in  an  overdense 
plasma,  which  has  a  higher  electron  density  ne  more  than  the  cutoff  density.  To 
eliminate  this  phenomenon  for  achievement  of  high-density  plasmas,  we  have  to  use 
external  magnetic  fields  to  modify  the  dispersion  relation  around  the  electron  plasma 
frequency  cope ,  at  which  ne  reaches  the  cutoff  density.  When  we  use  permittivity  s 


and  refractive  index  N  for  description  of  wave  propagation,  the  reason  why  the  wave 
cannot  propagate  in  an  overdense  plasma  is  that  N  is  imaginary  since  s  is  negative, 
via  the  equation  in  the  textbook  N  =  4s  ,  when  the  wave  angular  frequency  co  is  less 
than  &>pe . 

Such  a  negative  s  state  of  plasmas  is  useful  when  we  use  plasmas  for  control 
of  electromagnetic  waves.  Here,  the  waves  do  not  necessarily  propagate  inside  a  plasma. 
Several  authors  reported  effects  from  several  viewpoints  on  such  constituent  plasmas  as 
reflectors,4  absorbers,4’5  and  antennae6  of  electromagnetic  waves.  In  contrast  with  the 
other  materials  in  these  devices,  parameters  of  the  plasmas  are  fairly  flexible  to  be 
adjusted,  and  its  existence  itself  is  controllable.  Due  to  the  limitation  of  available  ne  in 
the  devices,  the  frequency  of  the  target  waves  is  up  to  THz  waves  so  far,  although  this 
limitation  can  be  removed  by  forth-coming  scientific  and  technological  progresses. 

Recently,  after  the  proposal  of  their  concepts,7  metamaterials8'9  have  attracted 
much  attentions  to  exhibit  extraordinary  features  such  as  negative  refractive  index,10 
perfect  lens  effects,11’12  and  cloaking.13  Here,  states  with  negative  s  are  quite 
important  to  make  N  negative,  via  the  formalism  used  in  the  metamaterial  area 
N  =  yfsyjji  and  with  a  negative-  /u  state,  where  //  is  (macroscopic)  permeability,  and 
are  inevitable  for  elaborate  control  of  N  for  cloaking.  In  an  adjacent  scientific  area, 
plasmonics  devices14  which  are  composed  of  metals  for  light  control  in  optics  also 
utilize  the  negative  s  states,  in  which  surface  plasmon  polaritons  are  excited  on  the 
metal  surfaces.  In  the  area  of  industrial  application  of  plasma  processing,  states  with 
negative  s  are  effective  for  large-area  plasma  generation  by  surface  waves,  leading  to 
relatively  high  ne  , 15,16  and  other  possibilities  can  be  expected  for  future  wave 
controllers. 

This  report  demonstrates  such  recent  progresses  of  wave  controllers  composed 
of  plasmas  and  other  materials,  stimulated  by  the  expanding  researches  of  metamaterials. 
However,  although  the  technique  to  realize  negative-  ju  states  is  introduced  from  the 
proposal  in  the  metamaterial  area,7  other  unique  characteristics  that  have  not  been 
observed  in  that  area  are  found  and  shown  here.  So  far,  after  the  previous  researches  of 
a  wave  controller  on  the  microstrip  line17’18  and  photonic-crystal-like  structures,17’19"24 
we  reported  effects  of  negative-  s  plasmas  in  theoretical  studies  on  localized 
surface-wave  modes  in  chain  structures  of  microplasmas,24"26  and  dynamic  negative-  N 
materials  composed  of  plasmas  and  metallic  resonators.27’28  The  entire  area  of  our 
researches  is  summarized  as  “plasma  metamaterials.”29 

In  this  report,  we  demonstrate  the  functional  composites  of  plasmas  and 
metamaterials,  and  the  focusing  point  is  verification  of  features  on  negative-  s  plasmas 
previously  predicted  by  our  studies.  The  distinguished  point  from  the  ordinary 


metamaterial  researches,  reported  here,  is  the  fact  of  flexible  and  complex  negative  s 
as  well  as  nonlinearity.30 


II.  Result  1:  Arbitrary  complex  functions  for  electromagnetic  waves  by  using  a 
plasma  array  reinforced  with  metallic  components  (ITEM  0001) 

States  with  negative  s ,  which  we  can  achieve  using  plasmas  and  can  be  reinforced  by 
metamaterials,  are  quite  useful  for  flexible  waveguides  for  microwaves.  Here  we 
confirm  effects  of  negative  s  in  composites  of  gaseous  plasmas  and  solid-state 
metamaterials,  whose  potential  roles  are  wave  controllers  which  can  regulate 
attenuation  and  phase  shift  independently. 

To  detect  a  rapid  change  of  complex-value  scattering  ( S )  parameters,  both 
waves  incident  into  and  transmitted  from  a  device  under  test  (DUT)  are  monitored  with 
good  time  response  faster  than  the  characteristic  time  of  its  dynamic  properties. 
Furthermore,  since  the  S  parameters  are  complex  variables,  a  detector  is  required  to 
extract  changes  of  both  their  amplitude  and  phase.  Here,  we  built  up  a  detecting  system 
using  microwave  devices  which  can  be  described  in  ramp  circuit  elements.  This  is 
similar  to  a  network  analyzer,  but  its  time  response  is  enough  to  trace  dynamic  changes 
of  plasma  parameters. 

Figure  1  shows  a  microwave  circuit  designed  for  this  purpose.  All  of  the 
devices  and  coaxial  cables  are  applicable  to  microwaves  up  to  18  GHz.  A  signal 
generator  (Agilent  Technology,  83624B)  provides  microwaves  at  frequencies  swept 
gradually,  although  the 
frequency  was  kept  constant  for 
time- varying  DUTs.  Incident 
microwaves  are  split  into  a 
probe  wave  and  a  reference 
wave,  where  the  probe  wave 
proceeds  to  a  DUT  and  the 
reference  wave  directly  goes  to 
the  S21  detector;  the  reference 
wave  gives  the  invariant 
measure  both  in  amplitude  and 
phase  components.  A  reflected 
wave  from  the  DUT  is  also 


Reflected  signal  (for  Sn) 


Fig.  1.  Microwave  circuit  used  for  dynamic  scattering 
parameters  Sn  and  S21. 


detected  at  Sn  detector  through  the  directional  coupler. 

To  detect  a  wave  signal  as  a  complex  value,  both  Sn  and  S21  signals  are 
filtered  in  a  hybrid  divider,  in  which  a  wave  is  split  into  two  waves  and  one  of  them 
suffers  phase  shift  of  90  degrees;  using  the  signals  in  phase  and  in  quadrature  phase,  we 
can  derive  both  amplitude  and  phase  components  of  the  wave.  We  refer  to  this  technique 
as  In-phase  Quadrature-phase  (IQ)  method  in  this  article.  We  note  that  a  conventional 
network  analyzer  has  a  processor  unit  that  derives  amplitude  and  phase  of  a  wave 
signals;39,40  the  unit  supports  precise  measurement  for  a  various  DUT,  but  requires  some 
processing  time  which  makes  it  difficult  to  catch  up  with  dynamic  change  of  DUT.  The 
experimental  results  shown  later  in  this  article  cannot  verify  its  applicability  for  every 
DUT,  although  it  works  well  for  our  DUTs,  i.e.,  with  -4  < Re(N)  <  4 ,  0<Im(N)<5, 
and  a  dynamic  property  with  the  characteristic  time  more  than  0.1  ps. 

As  usually  performed  in  a  network  analyzer,  a  calibration  method  to  minimize 
error  signals  is  required;  the  measurement  plane  at  each  port  of  a  two-port  device  is  not 
always  the  reference  plane  due  to  connection  with  series  and  parallel  floating  impedance. 
We  performed  this  procedure  according  to  the  literature.39,40  Furthermore,  another  error 
originates  from  multi  reflection  on  the  interface  planes  at  the  DUT  front  and  end,  which 
is  observed  as  periodic  humps  along  the  frequency  axis.  This  type  of  errors  can  be 
removed  through  series  Fourier  and  Inversed-Fourier  transform  procedure;  changing 
time  evolutions  to  the  frequency  domain  and  removing  fixed  error  spectra,  we  obtained 
pure  signals  coming  from  the  source.  The  entire  procedures  were  checked  out  using 
devices  and  waveguides  whose  properties  were  well  known,  and  we  could  confirm  that 
detected  signals  were  sufficient  for  physical  analysis  at  frequencies  ranging  from  2  to  7 
GHz. 

To  prepare  for  a  specific  DUT,  two  spiral  wires  covered  by  insulator  layers 
formed  a  double-helix  structure,  which  leads  to  series  inductance-capacitance  resonance, 
and  its  array  work  macroscopically  as  a  negative-permeability  material  if  it  has  an 
external  short-cut  path  and  its  loss  component  is  negligibly  small.27, 28  We  note  that  size 
of  its  outer  diameter  was  1.4  mm,  much  smaller  than  the  wavelength  of  the  microwaves 
at  7  GHz  (-43  mm  in  the  free  space).  The  two  insulated  wires  simultaneously  worked  as 
an  electrode  of  a  dielectric-barrier  discharge.  Since  they  were  installed  in  a  glass 
capillary  which  held  a  He  or  Ne  flow,  plasmas  at  elevated  gas  pressure  were  generated 
inside  and  outside  of  the  double-helix  structure.  These  plasmas  made  e  a  complex 
value  whose  real  part  could  be  negative  if  ne  was  so  high  that  cope  was  well  beyond 
co .  This  array  composed  of  15  units  with  its  length  of  30  mm  was  put  on  the  surface  of  a 
coplanar  waveguide  with  the  length  around  60  mm,  and  here  we  consider  this  coplanar 
waveguide  with  the  array  as  a  DUT.  Both  ends  of  the  coplanar  waveguide  were 


connected  to  the  microwave  circuit  shown  in  Fig.  1  via  coaxial  connectors.  Detailed 
installation  and  configuration  of  the  DUT  is  shown  in  Fig.  2.  Since  the  array  of  the 
metamaterial  was  coupled  with  electric  and  magnetic  fields  of  the  propagating 
microwaves  on  the  coplanar  waveguide,  the  electric  fields  are  approximately  parallel  to 
the  glass  capillaries.27’40  More  rigorously,  in  this  setup,  we  do  not  expect  propagation  of 
surface  waves  but  performances  as  macroscopic  media  emerging  in  composites  of 
plasmas  and  metal  structures.  The  discharge  voltage  applied  on  the  wired  electrodes  was 
operated  at  5  kHz,  and  so  this  component  which  might  affect  the  microwave  waveguide 
could  be  removed  by  bias-T  devices  set  on  both  sides  of  this  DUT. 

We  performed  experiments  of  microwave  transmission  through  this  DUT  using 
the  microwave  circuit  described  earlier.  A.  In  addition  to  the  simple  negative-  s  effects, 
by  varying  gas  conditions  and  electric  power  for  plasma  generation  to  control  Im(^) 
simultaneously  and  independently  with  Re(^),  we  will  see  the  entire  properties  of 
plasmas  with  negative-  //  metamaterials.  Note  that  //  determined  by  the  double  helix 
metal  wires  took  a  negative  value  just  above  magnetic  resonance  frequencies  and  a 
positive  value  outside  the  resonances; 28  one  of  the  magnetic  resonance  frequencies  was 
around  4.29  GHz.  We  here  note  smallness  of  the  value  of  |S21|  (—  0.1)  of  this  DUT  in  the 
following.  Since  we  set  the  working  point  with  negative  // ,  the  DUT  includes 
significant  loss  by  enhancement  of  Im(N)  ;  outside  magnetic  resonances, 
Im(N)/Re(N)  <  0.05,  whereas  Im(Nj  ~  Re( N)  without  plasmas  at  4.29  GHz  around 
which  Re(//)  is  negative,  mainly  due  to  large  Im(//).  One  might  suggest  that  large 
Im(//)  will  arise  from  extraordinary  radiation  losses,  but  we  did  not  see  such  signs  in 
measurements  using  a  scanned  micro  monopole  antenna  around  the  DUT  and  the 
waveguide;  radiation  signals  were  several  times  lower  when  we  installed  this 
metamaterial  structures  on  the  coplanar  waveguide  around  the  frequency  of  negative  // , 
and  we  conclude  that  such  a 
loss  took  place  inside  the 
metamaterial  structure.  This 
loss  is  inevitable  when  we  use 
inductive-capacitive 
resonances,9  although  lossless 
magnetic  resonances  using 
dielectric  resonators40  can  be 
applied  to  future  plasma 
metamaterials. 

Microscopic  profiles  of 
electric  and  magnetic  fields  are 


complicated  in  similar  manners  to  those  in  usual  metamaterial  experiments,9  and  we 
cannot  specify  detailed  propagating  modes  on  the  coplanar  waveguide  with  the  DUT 
using  a  simple  and  conventional  technological  term.  However,  since  the  fundamental 
mode  on  the  coplanar  waveguide  have  electric  fields  which  are  almost  parallel  to  the 
capillaries  and  perpendicular  to  the  spiral  wires,  the  electric  fields  mainly  suffer  plasma 
effects.  On  the  other  hand,  the  magnetic  fields  resonate  with  the  wires  at  4.29  GHz  via 
the  inductance-capacitance  resonance.  As  a  whole,  effective  macroscopic  parameters  of 
the  media  are  determined  through  such  plasma  metamaterial  effects,  and  can  be 
analyzed  by  S21  and  Sn  ,9'43 

Specifically,  we  used  He  and  Ne  gases  to  control  the  working  point 
significantly  on  the  Re^j-Im^)  complex  plane.  He  and  Ne  have  different  cross 
sections  of  elastic  collisions,44  by  changing  gas  pressure,  we  can  change  the  working 
point  of  s  on  the  plane.  Note  that  the  generated  plasmas  described  here  are  collisional 
to  see  the  effects  of  the  varying  working  points  on  the  complex  plane;  surface  wave 
propagation,  which  is  possible  in  collisionless  plasmas  as  shown  in  Section  III,  can  be 
ruled  out. 

In  the  case  of  He  at  atmospheric  pressure  (100  kPa)  near  the  resonance  (at 
4.297  GHz),  we  showed  the  results  of  time  evolutions  of  the  discharge  signals  in  Fig.  3. 
Using  the  signal  of  phase  and  amplitude  of  S21  in  this  time  evolution  and  extracting 
the  first  phase  of  the  discharge  voltage  in  Fig.  3,  we  can  make  a  new  drawing  of  the 
varying  working  point  of  S21  on  the  complex  plane  in  Fig.  4(a).  The  plasmas  were 
generated  in  a  scheme  of  dielectric  barrier  discharge.  From  Fig.  3,  we  recognized  that 
the  maximum  ne  was  observed  just 
after  the  rising  time  of  the  discharge 
voltage  (  t  ~  5  ps),  and  that  ne 
decreased  monotonously  until  the 
falling  time  of  the  discharge  voltage  ( t 
~  20  ps);  the  axis  of  t  can  be 
converted  to  ne  .  That  is,  such  new 
drawings  reveal  dynamic  and  tunable 
properties  of  this  metamaterial,  which 
are  expressed  by  the  relation  among 
three  variables:  |S21|cos(arg(S21))  , 

|S21|sin(arg(S21))  and  ne. 

In  Fig.  4(a),  the  trajectory  was 
along  one  straight  line  on  the  plane, 
which  indicates  attenuation  of  the  waves 


Fig.  3.  Time  evolutions  of  S21  and  discharge 
signals  in  DUT  at  100  kPa  of  He. 


took  place  without  phase  shift.  Also  shown  in  Fig.  7(a),  plasmas  are  equivalently 
attenuators  against  microwaves  in  almost  all  cases  except  those  shown  in  Fig.  8(b)  and 
(c);  electron  collisions  against  neutral  particles  leads  to  Im(f),  which  may  result  in 
Im(N) . 

Figure  4(b)  shows  varying  S21  on  the  complex  plane  in  the  case  of  He  at 
lower  pressure  (27  kPa)  near  the  resonance  (at  4.297  GHz).  This  case,  we  observed  a 
similar  trajectory  to  that  in  Fig.  4(a),  but  some  swings  due  to  phase  shift  was  slightly 
found.  In  contrast.  Fig.  4(c)  shows  varying  S21  on  the  complex  plane  in  the  case  of  Ne 
at  lower  pressure  (27  kPa)  near  the  resonance  (at  4.297  GHz).  We  observed  clear  phase 
shift  in  the  time  evolution,  which  means  a  significant  change  of  Re(N)  occurred. 
Furthermore,  if  // ,  which  is  determined  not  by  plasmas  but  by  the  double  helix  metal 
wires,  changes  to  the  positive  value,  the  result  is  totally  different;  Fig.  4(d)  shows  S21 
on  the  complex  plane  in  the  case  of  Ne  at  lower  pressure  (20  kPa)  outside  the  resonance 
(at  6.0  GHz),  and  we  can  see  only  attenuation  without  significant  phase  shift. 

These  properties  observed  in  the  experiments  include  the  DUT,  the  coplanar 


(a)  (c) 


Fig.  4.  Time  evolutions  of  S21  on  complex  plane,  (a)  Near  resonance  at  100  kPa  of  He.  (b)  Near 
resonance  at  27  kPa  of  He.  (c)  near  resonance  at  20  kPa  of  Ne.  (d)  Off  resonance  at  27  kPa  of  Ne. 


waveguide,  and  the  space  on  the  upper  side  of  the  DUT.  To  derive  macroscopic 
parameters  of  the  metamaterial  installed  on  the  upper  side  of  the  DUT,  we  adapted  one 
theoretical  model  developed  for  layered  permittivity  structure.29’45  In  Fig.  4(c),  more 
than  a  2 n  rad.  swing  of  the  trajectory  was  observed.  Taking  account  of  the  fact  that  the 
length  of  the  DUT  was  ~38  mm,  it  is  concluded  thatRe(N)  of  the  metamaterial,  the 
composite  of  the  plasmas  and  the  magnetic  resonators,  reached  -1.5 — 3  when  the 
maximum  phase  shift  took  place;  a  dynamic  negative-  N  material  was  successfully 
observed. 

Here,  we  mention  about  variation  of  such  macroscopic  parameters  depending 

no  no  nn  a/: 

on  gas  conditions  and  // .  Figure  5  displays  ne  dependence  of  s .  ’  ’  ’  As  the 
collisional  frequency  becomes  lower,  both  changes  in  Re(<s-)  and  Im(£-)  get  larger.  On 
the  other  hand,  higher  collisional  frequency  makes  the  ratio  Im(f)/Re(f)  larger. 
Similar  tendencies  are  found  for 
R e(js)  and  lm{js )  .  These  facts 
indicate  that  the  change  of  electron 
elastic  collisions  by  selection  of  gas 
species  and  pressure  leads  to  wave 
attenuation  and  phase  shift,  which  are 
controllable  independently.  Furthermore, 
due  to  the  rotation  of  the  working  points 
on  the  complex  4s  by  product  with 
■yfju  ,  variation  of  ne  makes  either 
change  of  Re(N)  or  Im(N);46  this  is 
one  of  the  unique  features  predicted  for 
plasma  metamaterials,  and  was  verified 
in  Fig.  4  (c)  and  (d).  Simultaneous 
realization  of  flexible  waveguides  and 
variable  attenuators  with  controllable 
phase  shifts  is  not  shown  in  this  report, 
but  it  will  be  possible  in  future  schemes 
of  localized  surface  wave  propagation 
on  negative  refractive  index  materials47 
composed  of  plasmas. 
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Fig.  5.  Complex  permittivity  s  as  a  function  of 
electron  density  nB.  (a)  He  at  100  kPa.  (b)  He  at  27 
kPa.  (c)  Ne  at  20  kPa. 


III.  Result  2:  Metamaterial-induced  transparency  of  plasma  layer  (ITEM  0001, 

0002) 

When  we  arrange  microplasmas  around  1  mm  along  one  line  on  one  plane,  microwaves 
with  its  frequency  slightly  less  than  a>pe  can  propagate  along  this  chain  strucrture, 24,25 
similar  to  localized  surface  plasmon  polaritons.14  Since  microwaves  form  standing 
waves  around  an  individual  microplasma  in  a  microscopic  point  of  view,  we  expect  that 
the  waves  can  propagate  along  bent  or  diverged  chains,  changing  the  direction  of  its 
macroscopic  wave  vector. 

A.  Numerical  analysis  of  propagating  localized  surface  waves  in  various  patterns  of 
plasmas 

When  we  arrange  microplasmas  around  1  mm  along  one  line  on  one  plane,  microwaves 
with  its  frequency  slightly  less  than  a>pe  can  propagate  along  this  chain  strucrture,24,25 
similar  to  localized  surface  plasmon  polaritons.14  Since  microwaves  form  standing 
waves  around  an  individual  microplasma  in  a  microscopic  point  of  view,  we  expect  that 
the  waves  can  propagate  along  bent  or  diverged  chains,  changing  the  direction  of  its 
macroscopic  wave  vector. 

First,  we  performed  numerical  analysis  to  confirm  these  predictions.  The 
method  we  used  was  different  from  the  one  for  derivation  of  dispersion  relation  in  an 
infinite  periodic  structure;  the  finite  difference  time  domain  (FDTD)  method  for 
dispersive  media  in  the  Drude  model  was  applied  to  study  this  phenomena,  since  it  can 
deal  with  finite  and  specific  structures.  Its  formalism  was  described  in  our  previous 
report,23  and  the  essential  point  is  that  we  introduced  one  extra  equation  for  balance  of 
current  densities  of  microwave  fluctuations  in  addition  to  the  Maxwell  equations.  We 
investigated  in-plane  electric-field  modes  which  have  electric  fields  parallel  to  the 
calculated  two-dimensional  plane  with  perpendicular  magnetic  fields.21,24,31  We  assumed 
that  microwaves  are  launched  by  a  micro  dipole  antenna,  which  is  given  by  a  FDTD  cell 
with  fixed  electric  field  in  the  direction  of  the  y  axis  near  the  plasma  column  on  the 
edge  and  clarifies  the  coupling  point  to  that  plasma;  port  analysis  for  scattering  matrix  is 
not  performed. 

Figure  6  shows  a  case  of  a  simple  straight  chain  of  microplasma  columns  with 
microwave  launching  points  and  calculated  electric  fields  expressed  by  their  amplitudes. 
We  assumed  radial  distribution  of  ne  in  the  usual  zero-order  Bessel  function  J0  in  the 
model,  where  it  is  monotonous  with  peaking  in  the  center  and  boundary  condition  of 
ne  =  0  on  the  edge.  The  exposed  electric  field  has  an  y  component  Ey  with  no  x 


component  Ex ,  and  this  Ey  works  for  a  boundary  condition  of  the  Maxwell  equations 
and  a  source  for  near  field  in  the  y  direction.  The  calculated  results  shown  in  Fig.  6(b) 
and  (c)  indicated  that  electric  fields  are  localized  on  the  edge  of  the  plasmas,  and 
propagate  along  the  chain  structure,  as  described  in  the  following.  Since  the  excitation 
point  is  very  close  to  the  first  microplasma  column,  the  excited  Ey  near  fields  induce 
surface  waves  on  the  round  surface  of  this  plasma  column,  where  the  surface  waves 
have  electric  field  components  parallel  to  and  perpendicular  to  a  round  interface 
between  a  positive  s  layer  (outside  region)  and  a  negative  s  layer  (inside  plasma 
region). 3,32,33  The  surface  waves  propagate  around  this  surface,  changing  its  propagation 
direction  and  its  field  components  as  mixtures  of  EXJ  and  Ex  that  are  symmetrical  in 
magnitude.  When  they  go  to  the  other  end  of  the  round  surface  of  the  first  plasma 
column,  another  microplasma  column  is  just  in  front  of  them;  most  of  the  wave  energy 
is  coupled  to  surface  waves  on  the  next  microplasma  column.  From  macroscopic  points 
of  view,  the  exposed  Ey  near  field  on  the  bottom  edge  of  the  column  chain  excites 
wave  propagation  along  the  chain  in  the  y  direction.  There  are  very  little  component 
of  propagating  (plane)  waves  excited  from  the  antenna  and  directly  coupled  to  the 
plasma  columns,  because  an  excited  propagating  wave  has  a  Poynting  vector  in  the  x 
direction.  The  maximum  value 


of  ne(- 1.5 xlO13  cm-3)  in  the 
center  of  the  individual  plasma 
corresponds  to  the  local 
electron  plasma  frequency  ~34 
GHz.  That  is,  the  excited  wave 
frequency  (13.1  GHz)  is  well 
below  the  cutoff  density  and  the 
wave  cannot  propagate  inside 
the  bulk  region.  The  value  of  s 
reaches  -20- -30  in  the  center, 
and  that  on  the  region  of 
localized  electric  fields  of  the 
surface  waves  is  0--1,  which 
are  suitable  for  surface  wave 
propagation  on  interface 
between  the  positive  s  layer 
(outside  region)  and  the 
negative  s  layer  (inside 
plasma  region). 
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Fig.  6.  (a)  Numerical  model  of  straight  plasma  chain  used 
here,  (b)  Ex  at  13.1  GHz  calculated  numerically,  (c)  Ey  at  13.1 
GHz  calculated  numerically. 
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Fig.  7.  (a)  Numerical  model  of  L-like-angled  plasma  chain  used  here,  (b)  Ex  at  13.1  GHz  calculated 
numerically,  (c)  Ev  at  13.1  GHz  calculated  numerically. 

Figure  7  shows  wave  propagation  along  an  L-shaped  bent  line  of  the  chains. 
The  wave  successfully  changes  its  propagation  direction  along  the  plasma  chain.  Such  a 
bent  propagation  path  is  possible  since  the  wave  propagates  microscopically  along  the 
edge  of  an  individual  plasma  and  partially  jumps  to  the  adjacent  plasma  through  a 
sneaking  component  normal  to  the  plasma  surface.  Note  that  the  wavelength  of  the 
surface  wave  is  much  shorter  than  propagating  waves  in  free  space  due  to  its  resonance 
property.  We  previously  reported  a  T-junction  of  similar  plasma  chain  structures,29  and 
these  facts  indicate  that  such  waveguides  can  be  designed  only  with  plasma  chain 
structures. 

B.  Experimental  observation  of  localized  surface  waves  and  their  propagation 

To  verify  the  prediction  of  the  numerical  results  shown  earlier,  we  performed 
experiments  in  the  similar  situation  of  the  numerical  model.  We  generated  columnar 
plasmas  with  sufficiently  long  length  (~  36  mm)  in  capillary- shaped  holes  in  a  dielectric 
cube  (teflon,  s  —2.1).  Using  Ar  gas  at  500  Pa,  we  generated  plasmas  by  pulsed  low 
frequency  waves  at  370  kHz  whose  amplitude  reached  around  0.8  kV.  The  gas  pressure 
we  set  assures  that  plasmas  generated  here  are  almost  collisionless,  which  means  that 
Re(^)  varies  with  almost  no  change  of  Imff) . 

Figure  8(a)  displays  visible  emissions  from  the  generated  plasmas,  viewed 
from  one  end  of  the  capillary;  one  end  was  a  mesh  electrode  to  pass  the  emission  light 
and  assure  the  flat  grounded  potential  profile.  The  plasmas  were  well  confined  in  each 
capillary  without  any  leakage,  and  formed  a  straight  line  composed  of  seven  ignited 
columns;  plasmas  were  not  generated  in  other  surrounding  capillary  holes.  Now,  two 


micro  monopole  antennae,  both  of  which  were  directed  in  the  y  axis  in  Fig.  1(a),  were 
set  at  both  ends,  and  one  was  for  the  launcher  of  microwaves  supplied  from  a  signal 
generator  (Agilent  Technology,  83624B),  and  the  other  was  the  receiver  with  a  rectifier 
to  monitor  the  amplitude  of  the  wave.  The  launcher  antenna  is  a  micro  monopole 
antenna  with  the  length  of  1  mm  on  a  microwave  coaxial  cable,  which  enables  us  to 


assure  similar  plasma-field  coupling  to  the 
emits  a  sufficient  component  of  electric 
near  fields  in  the  direction  of  y  axis. 
Figure  8(a)  also  shows  a  detected  signal 
profile  when  the  transmitted  signal  was 
fairly  strong  (at  1.51  GHz).  Its  profile  was 
quite  similar  to  the  electric  field 
distribution  at  the  end  of  the  chain 
calculated  numerically,  as  shown  in  Fig.  6, 
for  example.  If  the  receiver  signals  came 
from  the  propagating  plane  waves  in  the 
y  direction,  we  could  detect  almost  no 
signals  since  the  launcher  and  the  receiver 
antenna  mainly  had  a  Ey  component,  but 
it  was  not  the  case  at  this  frequency. 
Surface  waves  numerically  investigated  in 
Figs.  6  and  7  have  electric  fields  of  both 
x  and  y  components,3,32,33  and  the 
detected  Ey  electric  fields  are  consistent 
with  this  predicted  features. 

Figure  8(b)  shows  frequency 
spectra  of  the  transmitted  waves.  Such 
multiple  peaks  of  the  spectrum  in 
“enhancement  factor”,  which  is  defined  as 
the  ratio  of  signals  between  those  with  and 
without  plasmas,  were  quite  similar  to  the 
dispersion  relations  predicted 

previously.  ’  ’  That  is,  spectra  without 
plasmas  indicate  microwave  propagation 
through  the  dielectric  with  periodic 
capillary  tubes,  and  spectrum  variation 
found  in  those  with  plasmas  in  comparison 


numerical  analysis  shown  Fig.  6  and  7;  it 


Launcher 


20 


1.51GHz 

< 

p 

© 

©/ 

o 

\  © 

©o©©© 

''  GOoOOO 

6-4-2  0 

2  4  6 

Antenna  position  (mm) 


(a) 


(b) 


Fig.  8.  (a)  Spatial  profile  of  detected  electric 
field  in  straight  plasma  chain  configuration,  (b) 
Frequency  spectra  in  similar  configuration. 


with  the  case  without  plasmas  reflects  microwave  propagation  along  chain-like  plasma 
columns.  Among  the  observed  spectra,  we  stress  here  that  the  spectrum  around  1.5  GHz 
shows  a  very  clear  evidence  of  microwave  chain-like  propagation,  as  described  in  the 
following.  The  spectrum  around  1.5  GHz  in  the  case  with  plasma  generation  shows 
clear  enhancement  in  signals  of  transmitted  waves,  whereas  the  signal  is  quite  low  in  the 
case  without  plasma  generation,  although  all  the  displayed  signals  in  Fig.  3  are  more 
than  10  times  larger  than  the  noise  levels;  the  signals  below  the  maximum  noise  level  (< 
0.1  in  units  of  the  y  axis  in  Fig.  3)  are  omitted  from  the  plotted  data.  On  the  other 
hand,  other  two  peaks  in  the  case  with  plasma  generation  above  2  GHz  are  very  close  to 
the  peaking  signals  without  plasmas,  and  confirmation  of  chain-like  propagation  is 
uncertain.  This  is  partly  because  (plane)  waves  may  propagate  in  free  space  in  these 
higher  frequencies;  since  we  cannot  complete  ideal  settings  in  all  cases,  especially  in 
high  frequency  regions,  the  monopole  antenna  may  excite  waves  with  some  components 
of  propagation  in  the  y  direction,  and  such  wave  components  will  obscure  signals  of 
chain-like  propagation.  The  intensity  variation  in  the  frequency  spectra  above  2  GHz 
may  come  from  standing  waves  between  the  launcher  and  the  receiver  antennae.  On  the 
other  hand,  below  2  GHz,  the  signals  in  the  cases  without  plasmas  were  so  low,  which 
indicates  sufficient  elimination  of  wave  propagation  in  free  space.  In  fact,  in  this 
frequency  range,  half-wavelength  standing  waves  cannot  exist  due  to  short  distance 
between  the  antennae  (~  20  mm);  most  of  the  detected  signals  in  the  cases  without 
plasmas  are  expanding  near-field  ones.  These  facts  support  experimental  confirmation 
of  the  similar  waves  predicted  in  the  numerical  models  for  the  spectrum  around  1.5  GHz. 
Another  point  we  can  understand  from  this  data  is  that  ne  might  be  lower  than  the 
cases  investigated  numerically;  possible  values  of  ne  would  be  the  order  of  1012  cm-3 
with  the  maximum  local  values  of  ®  /In  -  10-20  GHz. 

The  enhancement  factor  shown  in  Fig.  8(b)  is  not  always  effective  on 
quantitative  estimation,  and  it  is  not  equal  to  physical  amplification  rate  due  to  the 
above  mentioned  influences  of  other  signal  mixing.  However,  from  this  factor,  we  can 
certainly  recognize  frequency  spectra  whose  widths  are  quite  small. 

Note  that  such  behaviors  of  surface  wave  propagation  on  negative-  s 
microplasma  chains  were  not  well  investigated  before  this  report.  The  surface  waves 
localized  on  individual  plasmas  and  connected  to  adjacent  plasma(s)  along  the  chain  are 
different  from  near  fields  which  are  simply  explained  by  the  Poisson’s  equation  and  no 
sensitive  dependence  on  wave  frequency.  Similar  wave  guiding  effects  have  been 
reported  as  surface  plasmon  propagation  on  nano-metallic  particle  chains14’  34'38,  but 
they  are  in  the  photon  range  and  with  different  spatial  s  profiles.  Here  we  show 
theoretical  predictions  on  details  of  the  surface  waves  on  plasma  chains  and  their 


experimental  evidences  on  localization  both  in  space  and  frequency  domain,  which  are 
only  explained  by  propagation  of  such  localized  surface  waves. 


IV.  Result  3:  Exploratory  application  to  high-power  microwave  control  by  plasma 
metamaterials  (ITEM  0002) 

When  energy  of  propagating  electromagnetic  waves  in  a  negative-permeability  space 
becomes  high,  the  waves  generate  plasmas,  and  the  plasmas  change  electric  permittivity 
in  the  corrresponding  space;  this  is  a  nonlinear  evolution,  and  should  be  treated  in  a 
self-consistent  analysis.  Here,  we  investigated  such  an  effect  theoretically  and 
experimentally,  and  they  generated  overdense  plasmas  successfully  and  part  of  the 
microwaves  can  transmit  in  a  nonlinear  scheme.  Theoretical  analysis  predicts  that 
high-density  plasmas  with  negative  permittivity  can  form  via  saddle-node  bifurcations 
within  an  adequate  electric  field.  To  confirm  theoretical  predictions,  using 
metamaterials  with  negative  permeability  achieved  by  magnetic  resonances,  we  inject 
microwaves  with  several  hundreds  of  watt  into  a  waveguide  filled  with  low-pressure 
discharge  gases,  and  we  confirmed  transmission  of  switchable  high-power  2.45  GHz  as 
well  as  overdense  plasmas  with  negative  permittivity. 

A.  Theoretical  analysis  of  nonlinear  system 

To  maintain  complete  self-consistency  between  the  microwave  energy  and  ne  in  a 
generated  plasma,  numerical  analysis  including  a  fluid  model  with  particle  balance  and 
Maxwell  equations  is  preferred,  but  it  is  rather  difficult  to  extract  underlying  physics 
since  calculated  results  obtained  in  numerical  analysis  is  fairly  complicated.  Here,  we 
use  the  method  reported  by  our  group  recently30,  reviewed  briefly  in  the  following.  The 
electric  fields  of  both  propagating  and  evanescent  waves  at  spatial  position  r  and  time 
t  have  phase  terms  including  N ,  expressed  using  wavenumber  k  and  frequency  of 
the  waves  co/ln  is  frequency  of  the  waves.  N  is  defined  in  regimes  of 
metamaterials  by  // ,  which  is  a  macroscopic  value  determined  by  synthesized  effects 
of  array  of  metamaterial  components,  and  s ,  which  is  a  microscopic  value  determined 
by  collective  motions  of  electrons  in  plasmas  and  given  in  the  Drude  model  using 
electron  plasma  frequency  cope ,  vacuum  permittivity  e0 ,  and  electron  mass  me .  We 
assume  collisionless  plasmas.  If  plasmas  are  generated  by  electric -field  drift  motions  of 
electrons  and  recombined  without  spatial  transport,  electron  continuum  equation  from 


the  particle  balance  expressed  with  electron  mobility  r/e  and  recombination  coefficient 
/?,  which  includes  the  terms  of  production  rate  and  recombination  rate.  Ionization 
coefficient  $¥alpha$  is  expressed  with  two  constants  C  =  29.2  and  D  =  26.6  for  Ag 
gas39  when  we  use  a ,  p  and  E  in  the  units  of  cm"1,  Torr  and  V/cm,  respectively.  In 
a  steady  state,  ne  is  as  a  function  of  E . 

This  system  includes  E ,  s  and  ne ,  and  we  solve  this  nonlinear  system  by 
comparing  two  plots  of  £  as  a  function  of  s  at  each  spatial  position.  We  use  similar 
numerical  results  reported  in  Ref.  30.  Then,  we  modify  the  parameter  of  Ar  gas  pressure 
to  100  Pa,  which  was  different  from  that  in  Ref.  30  and  similar  in  the  experiments 
shown  below,  and  calculate  £ .  Finally,  we  will  obtain  bifurcation  diagrams  which 
indicate  e  transitions  when  we  change  injection  power  of  microwaves  P{  or  electric 
field  at  the  source  point  E1 .  Note  that  ne  is  uniform  in  the  metamaterial  region  with 
slight  spatial  modulation  that  reflects  discrete  unit  structures  in  the  metamaterials,  and 
loss  mechanisms  of  electrons  are  recombination  in  a  given  position  at  which  they  are 
generated. 

First,  we  explain  the  model  and  the  results  of  theoretical  analysis  using  Fig.  9. 
Figure  9  (a)  shows  numerical  model  calculated  by  the  finite  difference  time  domain 
method.  From  the  source,  continuous  sinusoidal  waves  at  2.45  GHz  was  launched,  and 
the  metamaterial  region  with  150x200  mm2  had  a  specific  //  =  -1.  s  and  //  in 
the  surrounding  region  are  both  unity.  By  changing  s  in  the  metamaterial  region,  we 
can  derive  stable  values  of  the  local 
electric  field.  Note  that  near  fields 
around  the  source  point  induce  wave 
propagation  inside  the  metamaterial 
region  since  the  point  source  was  so 
close  to  the  region,  by  10-mm 
distance. 

Figures  9  (b)  and  (c) 
demonstrate  examples  of  numerical 
results.  When  ne  is  low  and  s  is 
positive  (Fig.  9  (b)),  N  is 
imaginary  since  //  is  negative.  As 
a  result,  the  waves  cannot  propagate 
inside  the  metamaterial  region.  On 
the  other  hand,  when  the  plasma  is 
overdense  and  s  is  negative  (Fig.  9 
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Fig.  9.  (a)  Numerical  model  used  here,  (b)  Electric 
fields  calculated  numerically  for  s  =  0.2.  (c)  Electric 
fields  calculated  numerically  for  s=  -6.0. 


(c)),  they  can  propagate  inside  the  region.  These  features  are  completely  on  the  contrary 
to  conventional  cases  with  positive  //  one  of  which  is  typically  shown  in  Fig.  9  (b). 
When  both  s  and  //  are  negative,  the  phase  velocity  is  reversed  although  the 
Poynting  vector  is  forward  from  the  wave  source  point  .29 

After  calculation  of  wave  propagation  in  various  cases  of  s ,  we  obtained  local 
electric  fields  En  at  several  points  along  the  propagation  path.  Figure  10  (a)  shows  the 
local  electric  fields  calculated  numerically.  Due  to  the  specific  geometrical  effects,  we 
recognize  data  scattering  irregularly,  but  roughly  we  can  estimate  dependence  of  the 
fields  by  solid  curves  as  a  function  of  e .  In  all  cases,  the  fields  are  almost  the  same 
when  e  is  fairly  negative.  On  the  other  hand,  at  the  points  apart  from  the  microwave 
source,  we  observe  almost  no  electric  fields  when  e  is  positive.  At  the  points  that  are 
very  close  to  the  microwave  source,  the  fields  are  strong  due  to  evanescent  waves  even 
in  the  cases  with  positive  s . 

Electric  fields  are  also  calculated  from  the  particle  balance  equation,  shown  as 
Er  in  Fig.  10  (b).  We  note  that  En  is  expressed  in  the  arbitrary  unit  and  linearly 
depends  on  E; .  That  is,  to  see  variation  of  ne  and  s  as  a  function  of  Ei ,  the 
solutions  are  given  by  crossing 
points  of  Er ,  which  are  the  fixed 
values,  and  En  with  varying  E1 . 

Figure  1 1  shows  s 
dependences  on  the  electric  field  at 
the  source  point  of  microwaves  E1 
which  is  directly  dependent  on  P{ . 

At  the  position  of  y  =  100  mm, 
evanescent  waves  with  near  fields 
are  dominant,  and  we  observe  only  a 
small  bifurcation  around  s  ~  0  and 
hysteresis  takes  place  with  a  small 
difference  of  E1  (0.3-0.4  in  the 
arbitrary  unit  in  Fig.  11).  At  the 
position  of  y  =  120  mm, 

propagation  waves  join  evanescent 
waves  which  are  depressed  as  the 
position  becomes  apart  from  the 
edge  of  the  metamaterial  region.  As 
a  result,  clearer  bifurcations  are 


(a) 


Fig.  10.  (a)  Electric  fields  calculated  numerically  on 
wave  propagation  for  various  e.  (b)  Electric  fields 
calculated  analytically  from  paraticle  balance  for 
various  f.. 


observed,  and  two  saddle-node 
bifurcation  points  at  E;  ~  1  and  0.6 
can  be  recognized.  That  is,  as  the 
working  point  goes  up  to  E;  ~  1, 
s  remains  around  +1.  Then,  at  the 
bifurcation  point  around  Ei  =1,  s 
jumps  up  to  — 20,  which  indicates 
high-  ne  plasma  generation  with  ne 
~  1012  cm"3.  After  the  high-  nc 
plasma  generation,  as  the  working 
point  goes  down  to  E;  ~  0.6,  the 
plasmas  are  in  overdense  states  with 
negative  s ,  and  then,  s  jumps  down  to  ~  +1  around  E{  =  0.6.  In  the  case  at  y  = 
200  mm,  around  which  the  propagating  waves  are  dominant,  significant  bifurcations 
with  larger  hysteresis  are  observed.  In  the  cases  at  y  =120  mm  and  200  mm,  no 
plasma  generation  is  expected  at  positive  e  and  the  plasmas  are  always  overdense. 

B.  Experimental  verification  of  enhanced  wave  transmission  and  overdense  plasma 
generation 

To  verify  the  above-mentioned  theoretical  predictions,  we  generated  plasmas  at  low 
pressure  using  metamaterials  with  macroscopic  negative  ju .  The  experimental  setup  is 
shown  in  Fig.  12.  The  conventional  rectangular  waveguide  at  2.45  GHz  was  installed  in 
a  vacuum  chamber,  which  was  filled  with  Ar  at  100  Pa.  This  gas  pressure  assures 
generation  of  collisionless  plasmas.  We  note  that  waveguide  structure  was  almost  kept 
both  outside  and  inside  the  vacuum  chamber,  only  with  small  discontinuity  at  the 
vacuum  window  made  of  quartz  glass  located  at  y  =0. 

We  set  an  array  of  double  split  ring  resonators7  as  a  negative-  //  metamaterial 
for  the  region  from  y  =  0  to  80  mm;  it  consisted  of  designed  copper  films  via  a  wet 
etching  process  from  plane  copper-covered  glass-epoxy  substrates.  By  derivation  of  the 
scattering  parameters  with  the  vector  network  analyzer  (Anritsu  Corp.,  MS2028B)  and 
using  the  parameter  retrieval  method,43  we  could  evaluate  the  value  of  macroscopic  // . 
In  the  case  in  Fig.  12(b),  the  magnetic  fields  $H$  penetrated  the  rings  and  we  expected 
achievement  of  negative  //  near  their  magnetic  resonance,  and  ju  was  derived  as 
-0.5-0.  lj  at  2.45  GHz.  On  the  other  hand,  in  the  case  in  Fig.  12(c),  the  magnetic  fields 
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Fig.  11.  Bifurcation  diagrams  between  s  and  initial 
electric  field  at  microwave  source. 
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Fig.  12.  (a)  Experimental  setup  for  high-power  microwave  injection  into  negative-//  metamaterial,  (b) 
Metamaterial  for  negative  ju.  (c)  Metamaterial  for  positive  ju. 
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Fig.  13.  (a)  Visible  image  of  plasma  emission  for  negative  ju.  (c)  Visible  image  of  plasma  emission  for 
positive  ju. 


H  never  penetrated  the  rings  and  we  could  not  expect  negative  //  since  there  were  no 
magnetic  resonances;  as  a  result,  //  was  derived  as  ~  +1  at  2.45  GHz. 

Then,  we  injected  2.45-GHz  microwaves  with  360  W  as  a  forward  power  with 
25-psec  width  from  a  high-power  microwave  amplifier  (Kyoto-Micro-Densi, 
MA-02400C)  into  the  vacuum  chamber.  Figure  13(a)  shows  visible  emission  of  plasmas 
generated  in  the  metamaterial  region  of  Fig.  12(b).  The  images  seemed  to  be 
inhomogeneous  due  to  the  view  angle,  but  they  were  uniform  on  the  cross  section  of  the 
metamaterial.  From  emissions  near  the  right  edge,  we  noticed  that  intensified  emissions 
were  near  the  entrance  of  the  metamaterial  region,  and  they  decreases  gradually  as  the 
microwaves  propagated.  In  comparison,  when  we  use  another  metamaterials  structure, 
shown  in  Fig.  12(c),  whose  DSRRs  are  perpendicular  to  those  in  Fig.  12(b),  we 


observed  very  weak  and  unstable  emissions 
as  shown  in  Fig.  13(b).  We  note  that  we 
detect  no  signals  above  the  detection  limit 
of  single  Langmuir  probes  in  the  case 
shown  in  Fig.  12(c). 

Now  we  are  concentrated  on 
results  of  the  negative  p  metamaterial 
shown  in  Fig.  12(b).  Figure  14(a)  shows 
time  evolutions  of  transmitted  power 
detected  by  the  directional  coupler  on  the 
other  side  of  the  entrance.  Without  plasma 
generation,  we  did  not  detect  any  signals, 
unlike  the  forward-  and  the  reflected-power 
signals.  However,  when  the  plasmas  were 
generated,  the  transmitted  power  was 
detected,  as  shown  in  Fig.  14(a),  and 
increased;  this  is  the  clear  evidence  for 
switchable  passive  microwave  transmitters 
based  on  this  plasma  metamaterial 
configuration.  Such  phenomena  were  not 
the  case  of  microwave  plasma  generation  in 
conventional  schemes;  usually,  after 
plasma  generation,  the  transmitted  power 
decreases  due  to  power  dissipation.  In  our 
case,  the  change  of  N  from  the  imaginary 
state  with  positive  s  and  negative  //  to 
the  real  state  with  a  double  negative  system 
induced  enhancement  of  transmittance. 

Figure  14  also  shows  ion 
saturation  current  7is  at  bias  voltage  of 

-20  V  measured  by  a  conventional  single 
probe  at  y  =10  mm;  this  signal  coincided 

with  the  power  signals.  By  changing  bias 
voltage,  we  obtained  probe  current-voltage 
curves  in  a  few  tens  of  shots.  Figure  14(c) 
shows  evolution  of  ne  and  electron 
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Fig.  14.  Time  evolutions  of  (a)  microwave  signals, 
(b)  probe  signal,  and  (c)  plasma  parameters. 
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Fig.  15.  Electron  density  ne  as  a  function  of 
forward  microwave  power.. 


Electron  temparature  T  (eV)  (Rear  directional  coupler)  (arb.  unit) 


temperature  Te .  The  solid  line  indicates  the  well-known  cutoff  density  ( 7.4  x  1010  cm'3) 

for  waves  in  the  TEM  mode  or  propagating  in  an  infinite  plasma,  and  the  dashed  line 

indicates  the  cutoff  density  (5.1x  1010cm'3)  for  those  in  the  TEi0  mode  or  propagating 
in  an  plasma  occupying  the  rectangular  waveguide.  Monitored  ne  is  always  well 

beyond  the  both  cutoff  densities;  overdense  plasmas  with  negative  s  were 
successfully  generated. 

Figure  15  shows  variation  of  ne  as  a  function  of  the  forward  power  Pf . 
When  Pf  was  quite  low  (less  than  40  W),  we  did  not  recognize  any  plasma  generation, 
and  when  it  was  higher  than  100  W,  we  always  observed  overdense  plasmas.  When  P( 
is  around  50  W,  we  detected  plasmas  in  some  discharge  shots  and  no  plasmas  in  the 
other  shots,  and  in  the  cases  of  plasma  generation,  we  always  observed  overdense 
plasmas;  we  note  that  s  =  0  at  the  cutoff  condition,  which  is  in  the  TEio  mode  since 
we  generated  plasmas  inside  the  waveguide.  This  experimental  results  indicate  two 
possible  states,  similar  to  bifurcated  solutions  of  ne  or  s  for  E;  =  0.6- 1.0  with  40 

mm  shown  in  Fig.  11. 
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